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Abstract. The modification of hard jets in dense matter has so far been described by four different for- 
malisms based on perturbative QCD (pQCD). In these proceedings, we compare the various approximations 
made in these different schemes, especially those regarding the structure of the medium through which 
jets propagate. Following this, we highlight some of the major differences in the various physical processes 
contained in the different approaches. 
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1 Introduction 

One of the major discoveries of the heavy-ion program 
at the Relativistic Heavy-Ion Collider (RHIC) has been 
the observed suppression of high transverse momentum 
(high px) hadrons when compared to the yield of simi- 
lar hadrons in p-p collisions (scaled up by the expected 
number of binary collisions) [HE]- In p-p collisions, such 
hadrons are formed in the fragmentation of high px jets 
produced in hard scatterings. The presence of a dense 
medium influences the space-time development of the par- 
tonic shower from such jets and in turn leads to a medium 
modification of the final fragmentation to hadrons [3l|4] . 

The presence of a hard jet introduces a large energy 
scale within the process and allows for a calculation of 
. the modification using the methods of perturbative QCD 
(pQCD). Following the early attempts of Baier-Dokshitzer- 
Mueller-Peigne-Schiff and Zakharov (BDMPS-Z) [HHZ1IE1 
[fjlfTOlITT] , such calculations have grown in both sophistica- 
tion and in the number of different observables that they 
are applied to. The majority of current approaches to the 
energy loss of light partons may be divided into four ma- 
jor schemes often referred to by the names of the original 
authors: 

— Higher Twist scheme (HT) [HUIMI11IIMIMI2] 

— Path integral approach to the opacity expansion by 
Armesto, Salgado and Wiedemann, (BDMPS-Z/ASW) 
[I*prjpTp2pgl|5^^ 

— Finite temperature field theory approach by Arnold, 
Moore and Yaffe (AMY) [24lf2^[251[27l[28] 

— Reaction Operator approach to the opacity expansion 
by Gyulassy, Levai and Vitev, (GLV) [29ll30p 3i p32(l33] 



a Present address: Department of Physics, Ohio State Uni- 
versity, Columbus, Ohio 43210, USA. 



All these schemes utilize slightly different approximations 
regarding the various scales involved in the calculation and 
somewhat different quantitative pictures of the medium. 

It will be demonstrated in the companion publication 
of Ref. [Mj, that, using these different formalisms to com- 
pute the medium modification of hard jets in an identical 
medium leads to rather similar predictions for experimen- 
tal observables. In these proceedings, we outline the vari- 
ous differences between the theoretical formulations of the 
different schemes themselves. In Sect. [21 we present a brief 
introduction to the basic formalism and the definition of 
a medium modified fragmentation function. In Sect. El we 
present a brief review of how the different schemes com- 
pute the single gluon emission spectrum. In Sect. [H we 
review how a single gluon emission spectrum is iterated 
in the different formalisms. We present concluding discus- 
sions in Sect. [5] 

2 Hard scattering and the medium modified 
fragmentation function 

In the collision of two heavy-ions, there occasionally oc- 
curs a hard scattering between two initial partons which 
leads to two back-to-back out-going partons with large 
transverse momentum. These encounter multiple scatter- 
ing in the produced medium leading to a modification 
of the final distribution of hadrons emanating from these 
partons. In the computation of this modified distribution, 
all schemes utilize a factorized approach where the final 
cross section to produce a hadron h with transverse mo- 
mentum pt (rapidity between y and y + dy) may be ex- 
pressed as a convolution of initial nuclear structure func- 
tions [Gq (x a ), Gfr(xb), initial state nuclear effects such 
as shadowing and Cronin effect are understood to be in- 
cluded] to produce partons with momentum fractions x a , xj,, 
a hard partonic cross section to produce a high transverse 
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momentum parton c with a transverse momentum p and 
a medium modified fragmentation function for the final 
hadron [D%(z)}, 



d 2 a h 1 

= - / dx a dx b G^(x a )G^{x b ) 



dyd 2 pT 7r 

d<Tab 



*cX 



D h c {z) 



dt 



(1) 



In the vicinity of mid-rapidity, z — pt/p and t = (p — 
x a P) 2 (P is the average incoming momentum of a nu- 
cleoli in nucleus A). The entire effect of energy loss is 
concentrated in the calculation of the modification to the 
fragmentation function. The four models of energy loss 
are in a sense four schemes to estimate this quantity from 
perturbative QCD calculations. 

While the terminology (medium modification) used to 
describe the change in the fragmentation function seems 
to indicate that the medium has influenced the actual pro- 
cess of the formation of the final hadrons from the partonic 
cloud, this is not the case. All computations simply de- 
scribe the change in the gluon radiation spectrum from a 
hard parton due to the presence of the medium. The final 
hadronization of the hard parton is always assumed to oc- 
cur in the vacuum after the parton, with degraded energy, 
has escaped from the medium. Note that some of the hard 
gluons radiated from the hard parton will also encounter 
similar "modification" in the medium and may endure vac- 
uum hadronization after escaping from the medium. Dif- 
ferences between formalisms also arise in the inclusion of 
hadrons from the fragmentation of such sub-leading glu- 
ons: whereas in approaches which compute the change in 
the distribution of final partons (such as AMY) or the 
change in the distribution of final hadrons (such as HT) , 
hadrons from sub-leading gluons are implicitly included, 
formalisms which compute the energy loss of the leading 
parton (such as ASW), do not include such sub- leading 
corrections. 

To better appreciate the approximation schemes, one 
may introduce a set of scales (see Fig. [T]): E or p + , the 
forward energy of the jet; Q 2 , the virtuality of the initial 
jet-parton; /i, the momentum scale of the medium and 
L, its spatial extent. Most of the differences between the 
various schemes may be reduced to the different relations 
between these various scales assumed by each scheme as 
well as by how each scheme treats or approximates the 
structure of the medium. In all schemes, the forward en- 
ergy of the jet far exceeds the medium scale, E 3> \i. 



3 Single gluon emission and scattering in the 
medium 

The first step in energy loss calculations is to compute 
the effect of a single gluon emission off a hard jet in the 
medium. The major theoretical differences between the 
various schemes arise in this calculation. It is in this step 
that differing assumptions regarding the medium (in dif- 
ferent formalisms) are introduced. In the next section, the 
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Fig. 1. A schematic picture of the various scales involved in 
the modification of jets in dense matter. 



single emission kernel will be repeated to compute the ef- 
fect of multiple emissions. In most cases, this involves a 
certain phenomenological picture and also introduces fur- 
ther differences between the different approaches. 



3.1 Higher twist approach 

The origin of the higher twist (HT) approximation scheme 
lies in the calculations of medium enhanced higher twist 
corrections to the total cross section in Deep- Inelastic Scat- 
tering (DIS) off large nuclei [35 . In those calculations, the 
authors computed a certain class of power corrections to 
the total leading twist cross sections, which, though sup- 
pressed by powers of the hard scale Q 2 , are enhanced by 
the extent of the medium. In the case of high px hadron 
production one identifies and resums corrections to the 
single hadron inclusive cross section. 

One presupposes that the produced jet has a very large 
forward energy E which is much larger than its virtu- 
ality Q (which limits the transverse momentum of the 
radiated gluon, fcj_), which in turn is much larger than 
the characteristic momentum scale in the medium /x, i.e., 
£ > fci > /J. This hierarchy is then applied to the com- 
putation of multiple Feynman diagrams such as the one 
in Fig. [21 This diagram represents the process of a hard 
virtual quark produced in a hard collision, which then ra- 
diates a gluon and then scatters off a soft medium gluon 
with transverse momentum q± ~ \i prior to exiting the 
medium and fragmenting into hadrons. Even at the order 
considered, there exist various other contributions which 
involve scattering of the initial quark, off the soft gluon 
field, prior to radiation as well as scattering of the ra- 
diated gluon itself. All such contributions are combined 
coherently to calculate the modification to the fragmenta- 
tion function directly. 

The hierarchy of scales allows one to use the collinear 
approximation to factorize the fragmentation function and 
its modification from the hard scattering cross section. 
Thus, even though such a modified fragmentation function 
is derived in DIS, it may be generalized to the kinematics 
of a heavy-ion collision. Diagrams where the outgoing par- 
ton scatters off the medium gluons, such as those in Fig.[3J 
produce a medium dependent additive contribution to the 
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vacuum fragmentation function, which may be expressed 
as, 



AD^z.Q 2 ) = f 
Jo 



Q2 dk 



k\ 2vr 



(2) 



In the above equation, AP^j represents the medium mod- 
ified splitting function of parton i into j where a momen- 
tum fraction x is left in parton j. The new momentum 
fraction xl = k\j (2P~ p + x(\ — where the radiated 

gluon or quark carries away a transverse momentum k±, 
P~ is the incoming momentum of a nucleon in the nu- 
cleus and p is the momentum of the virtual photon. The 
medium modified splitting functions may be expressed as 
a product of the vacuum splitting function Pi—*? and a 
medium dependent factor, 



AR. 



P, 



(N? l)q 



2irC R (kl 



(Ql)) 



f(C,x L ). (3) 



Where, Cr is the representation dependent Casimir and 
N c is the number of colours. The mean transverse momen- 
tum of the soft gluons is represented by the factor (q\). 
The distance £ is the distance between the origin of the jet 
and the location of its scattering, which is limited by the 
length of the medium L. The function /(O^l) depends 
on the number of scatterings per radiated gluon included 
and encodes the in-medium interference effects such as the 
Landau-Pomeranchuck-Migdal effect [361137] . 

The factor q encodes the soft gluon field in the medium, 
off which the jet encounters multiple scattering. It is given 
as M 
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(4) 



x (F->(( + f/2, £J2)F°-(C - C+/2, -Cj./2)). 

The transport coefficient is normalized by fitting to one 
data point and a model such as a Woods-Saxon distribu- 
tion for cold matter or 3-D hydrodynamical evolution for 
hot nuclear matter is invoked for its variation with space- 
time location. The expectation ( ) is meant to be taken in 
the medium under consideration. Any space time depen- 
dence is essentially included in the implied expectation. 

The gluons which contribute to q do not have to be 
the entropy carriers of the system. In applications to cold 
nuclear matter, these gluons constitute the virtual gluon 
cloud inside the nucleons. In the case of a deconfined 
quark-gluon plasma, these may be the entropy carrying 



1 Throughout the HT portion of these proceedings, four- 
vectors will often be referred to using the light cone convention 
where x + = x° + x 3 and x~ = (x° - x 3 )/2. 




Fig. 2. A typical higher twist contribution to the modification 
of the fragmentation function in medium. 



gluons or virtual excitations within these degrees of free- 
dom. Which gluons the jet scatters off depends on the scale 
of the hard jet. It is immediately obvious from Eq. ([4} that 
q is a function of the jet energy p + . Note that p + is not 
integrated out. The actual dependence on p + depends on 
the medium in question. In the case of confined nuclear 
media, or a quark gluon plasma, the dependence is log- 
arithmic. There is also a logarithmic dependence on the 
virtuality of the jet which sets in due to radiative correc- 
tions to the definition in Eq. ((3]). Also, as demonstrated 
in Ref. [18] . q may even possess a tensorial structure if 
the medium is not isotropic. In the calculations of the 
current manuscript, both the dependence on the energy 
and virtuality of the jet will be ignored. The medium will 
be assumed to be isotropic. The values of q quoted should 
thus be considered as approximations to the full functional 
form. 



3.2 Opacity expansion approach 

Unlike the higher twist scheme, which is set up to directly 
calculate the final distribution of hadrons, opacity expan- 
sion approaches such as the Gyulassy-Levai-Vitev (GLV) 
scheme [29,30,31 , and the Armcsto-Salgado- Wiedemann 
(ASW) scheme p^l[20ll21U22l[23] were constructed primar- 
ily to deal with the problem of energy loss of the leading 
parton in dense deconfined matter. Both these schemes 
assume that the medium is composed of heavy almost 
static color scattering centers which are well separated, 
in the sense that the mean free path of a jet A the 
color screening length of the medium [4]. The opacity of 
the medium n, which constitutes the expansion parameter 
of these calculations, quantifies the number of scattering 
centers seen by a jet as it passes through the medium, 
i.e., n — L/X, where L is the thickness of the medium. 
The difference between the two approaches of the GLV 
and the ASW arise from how these tend to expand in n. 
In the GLV formalism, one constructs a recursive opera- 
tor expansion in opacity, whereas in the ASW approach 
a path integral over opacity is formulated. The solution 
of the recursive operator approach in the GLV allows for 
an order-by-order expansion in opacity. The path-integral 
in the ASW approach has been solved analytically in two 
limits: the one-scattering approximation, equivalent to a 
first order in opacity calculation in the GLV approach and 
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in the multiple scattering approximation where all orders 
in opacity have been resummed. In this article (as well as 
in the companion [34] where results of calculation will be 
compared with experimental data), the focus will lie on 
the path integral approach of ASW. 

The path integral approach for the energy loss of a 
hard jet propagating in a colored medium was first intro- 
duced in Ref. 8 . It was later demonstrated to be equiv- 
alent to the well known BDMPS approach [SHUE] m the 
multiple scattering limit. ASW represents the current, most 
widespread, variant of this approach. In this scheme, a 
hard, almost on-shell parton traversing a dense medium 
full of heavy scattering centers will engender multiple trans- 
verse scatterings of order [i -C p + . It will in the process 
split into an outgoing parton and a radiated gluon which 
will also scatter multiply in the medium. The radiated 
gluon, induced by the multiple scattering, has a transverse 
momentum k± > fi (different from the HT approach) . The 
propagation of the incoming (outgoing) partons as well as 
that of the radiated gluon in this background color field 
may be expressed in terms of effective Green's functions 
[G(r±, z; r± ', z') (for quark or gluon)] which obey the ob- 
vious Dyson- Schwinger equation, 

G(r±,z; r'x, z') = G (r±,z; r'±, z') 
- i Jd(Jd 2 xG (r ± ,z; x, ()A (x, ()G(x, C; r'±, z'),(5) 

where, Go is the free Green's function and Aq represents 
the color potential of a scattering center in the medium. 
The solution for the above interacting Green's function in- 
volves a path ordered Wilson line which follows the poten- 
tial from the location [r±(z'),z'] to [ri(z),z]. Expanding 
the expression for the radiation cross section to order A^ n 
corresponds to an expansion up to n th order in opacity. 

Taking the high energy limit and the soft radiation ap- 
proximation (x <C 1), one focuses on isolating the leading 
behavior in x that arises from the large number of interfer- 
ence diagrams at a given order of opacity. As a result of the 
approximations made, one recovers the BDMPS condition 
that the leading behavior in x is contained solely in gluon 
re-scattering diagrams. This results in the expression for 
the inclusive energy distribution for gluon radiation off an 
in-medium produced parton as [20] . 



conjugate whereas yi and yi represent the longitudinal po- 
sitions. The density of scatterers in the medium at location 
£ is n(() and the scattering cross section is a(r). In this 
form, the opacity is obtained as J n(£)<i£ over the extent 
of the medium. 

Numerical implementations of this scheme have fo- 
cused on two separate regimes. In one case, cr(r) is re- 
placed with a dipole form Cr 2 and one solves the harmonic 
oscillator like path integral. This corresponds to the case 
of multiple soft scatterings of the hard probe. In the limit 
of a static medium with a very large length, one obtains 
the simple form for the radiation distribution 38J, 

dl_ 2a s C R { \f^> 2 for w < ^c, 
"fa" n fora;>u, c . 

Where lu c = /rfCC<?(C) is called the characteristic fre- 
quency of the radiation. Up to constant factors, this is 
equal to mean energy lost in the medium ((E)) i.e., lo c — 
2(E) /(a s Cn). For a static medium, the integral defining 
lo c may be performed to obtain uj c = qL 2 /2, where L is 
the length of the medium and q is the jet transport coef- 
ficient, defined as the transverse momentum picked up by 
a hard jet per unit length. In actual numerical implemen- 
tations, the mean q (or the q at a well defined location 
and time) is the only tunable parameter when comparing 
with experimental data. For a dynamical medium of finite 
extent, the characteristic frequency and the overall mean 
transverse momentum gained by the jet (qL) will have to 
be estimated based on an Ansatz for the space time dis- 
tribution of the transport parameter q (see Ref. |34] for 
further details). 

In the other extreme, one expands the exponent as a 
series in na; keeping only the leading order term corre- 
sponds to the picture of gluon radiation associated with 
a single scattering. In this second form, the analytical re- 
sults of the ASW scheme formally approach those of the 
GLV reaction operator expansion [35]. In either case, the 
gluon emission intensity distribution has been found to 
be rather similar, once scaled with the characteristic fre- 
quency in each case. 



3.3 Finite temperature field theory approach 
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(6) 



where, as always, k± is the transverse momentum of the 
radiated gluon and xp + is its forward momentum. The 
vectors y and u represent the transverse locations of the 
emission of the gluon in the amplitude and the complex 



In this scheme, often referred to as the Arnold-Moore- 
Yaffe (AMY) approach, the energy loss of hard jets is con- 
sidered in an extended medium in equilibrium at asymp- 
totically high temperature T — > oo. Due to asymptotic 
freedom, the coupling constant g — > at such high tem- 
peratures and a power counting scheme emerges from the 
ability to identify a hierarchy of parametrically separated 
scales T ^> gT ^> g 2 T etc. In this limit, it becomes pos- 
sible to construct an effective field theory of soft modes, 
i.e., p ~ gT by summing contributions from hard loops 
with p ~ T, into effective propagators and vertices [39j . 

One assumes a hard on-shell parton, with energy sev- 
eral times that of the temperature, traversing such a medium, 
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Fig. 3. A typical cut diagram in the AMY formalism. 



undergoing soft scatterings with momentum transfers ~ 
gT off other hard partons in the medium. Such soft scat- 
terings induce collinear radiation from the parton, with a 
transverse momentum of the order of gT. The formation 
time for such collinear radiation ~ l/(g 2 T) is of the same 
order of magnitude as the mean free time between soft 
scatterings [35]- As a result, multiple scatterings of the 
incoming (outgoing) parton and the radiated gluon need 
to be considered to get the leading order gluon radiation 
rate. One essentially calculates the imaginary parts of in- 
finite order ladder diagrams such as those shown in Fig. [3l 
this is done by means of integral equations [26] . 

The imaginary parts of such ladder diagrams yield 
the 1 — > 2 decay rates of a hard parton (a) into a ra- 
diated gluon and another parton (6) Tg . These decay 
rates are then used to evolve hard quark and gluon dis- 
tributions from the initial hard collisions, when they are 
formed, to the time when they exit the medium, by means 
of a Fokker-Planck like equation [27], which is written 
schematically as, 



dPgjp) 

dt 



dk 



b,c 



Pa(p) 



dr b a c ( P ,k) 

dkdt 



k) 



dr b ac { P + k lP ) 

dkdt 



(8) 



The use of an effective theory for the description of 
the medium and the propagation of the jet, makes this 
approach considerably more systematic than the two pre- 
vious approaches: both the properties of the jet and the 
medium are described using the same hierarchy of scales. 
It remains the only approach to date which naturally in- 
cludes partonic feedback from the medium, i.e., processes 
where a thermal quark or gluon may be absorbed by the 
hard jet 0. In contrast to ASW and HT, this approach 
also (naturally) includes flavor changing interactions in 
the medium. Elastic energy loss may also be incorporated 
within the same basic formalism [41 j . However, since the 
AMY scheme assumes a thermalized partonic medium, its 
applicability is somewhat limited: It cannot compute the 
quenching of jets in the confined sector. As a result, en- 
ergy loss in cold confined nuclear matter as well as in 



2 While an attempt to include such effects in the higher twist 
formalism have been made in Ref. [40], these remain as phe- 
nomenological extensions and have not been included in this 
manuscript. 



the hadronic phase of heavy-ion collisions cannot be com- 
puted in this model. The off-shellness or virtuality of all 
jets is considered to be similar to that of hard partons in 
the medium, as a result, interference between vacuum and 
medium induced radiations is also not considered. 

In realistic calculations, the temperature of the medium 
is usually set by the underlying hydrodynamic simulation 
(see Ref. [31] for details). While in the HT or ASW for- 
malisms, an Ansatz is made for the one tunable parameter 
q and its relation to T, in the AMY formalism, q may be 
calculated directly from a knowledge of the temperature 
and the strong coupling constant g (or a s ). This is due 
to the precise picture of the medium used: that of a hot 
plasma of quarks and gluons. In realistic simulations, the 
coupling is, in principle, unknown and becomes the pri- 
mary fit parameter. This is then fit by comparison to one 
data point. 



4 Multiple gluon emissions 

In the preceding section, the effect of a single gluon emis- 
sion, stimulated by scattering in the medium, was consid- 
ered. In order to compute the final spectrum of hadrons, 
this single emission kernel has to be be repeated to ac- 
count for multiple gluon emissions and folded with a non- 
perturbative fragmentation function. Even in this proce- 
dure, the different schemes employ different methods: The 
HT scheme starts with a fragmentation function at a lower 
scale [i and evolves this distribution up to a higher scale. 
The ASW scheme, considers a finite energy lost by the 
leading parton in multiple unrelated events by means of 
a Poisson distribution and folds the outgoing parton with 
a vacuum fragmentation function with a shifted momen- 
tum fraction z. The AMY formalism considers the evolu- 
tion of an initial distribution of hard partons with time in 
the medium using the Fokker-Planck equation afforded by 
Eq.[S]and final also uses a vacuum fragmentation function 
with a shifted momentum fraction z. 



4.1 Higher twist scheme 



In subsection l3.H the medium modified fragmentation func- 
tion calculated in Eq. ([3]) included only one gluon emission 
in the medium. Any remaining gluon emissions occurred 
in the vacuum and were included in the renormalization of 
the vacuum fragmentation function. Unlike the remaining 
formalisms, the results from just the single gluon emission 
in the medium yield a medium modified fragmentation 
function and are already comparable with experiment. 

In reality, one expects multiple emissions to occur in 
the medium, followed by escape into the vacuum and fur- 
ther emissions in the vacuum. One starts with a vacuum 
fragmentation function at a low scale fJ.i ow and insists 
that the parton exits the medium with a certain virtuality 
/z. Emissions from the scale \x\ ow up to the scale /j may 
be included by using the standard vacuum Dokshitzcr- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equa- 
tions gaiiiiiB!. 
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Emissions in the medium account for the remaining 
evolution from the scale \x up to the scale Q. To compute 
this in-medium evolution, the medium modified fragmen- 
tation function from single gluon emission in Eq. ([3]) is 
now generalized to an evolution equation in virtuality of 
the propagating parton (see Ref. [45] for details), i.e, 



dD%(z,M 2 ,p+)f c f i 
dlog(M 2 ) 



i O 

■= / dy J dQPi >r<w , 

Ci 



(N 2 C - l)g(C) 
■2nC R (ki + (qi)) 



x f{C,x Ll y)D h q ( ^,M 2 ,q-y 



it 



(9) 



The initial conditions to this differential equation are 
provided by the fragmentation functions at the scale \x. 
The final resulting medium modified fragmentation func- 
tions includes both vacuum and in-medium induced emis- 
sions from the scale Q down to the scale fx. Further emis- 
sions occur solely in the vacuum. This medium modified 
fragmentation function may now be convoluted with the 
cross section to produce a hard parton [as in Eq. ([1])] 
to find the final distribution of hadrons. Computation of 
the medium modified fragmentation function in an evolv- 
ing medium such as the deconfined matter formed in a 
quark-gluon plasma involves further calculational details 
presented in Ref. 34j. As the HT formalism is setup to 
directly calculate the final modified fragmentation func- 
tion, it offers the simplest and most direct extension to 
the study of multi-hadron observables [46] , 



4.2 Opacity expansion scheme 



In subsection 13. 2\ the opacity expansion was used to cal- 
culate the differential spectrum for single gluon radiation 
from a hard parton. The calculations were carried out in 
the soft gluon limit i.e., lu — > 0. The next step is to calcu- 
late the probability for the leading hard parton to radiate 
off a finite energy AE = eP + . After losing this energy, the 
degraded hard parton escapes the medium and fragments 
in vacuum into a shower of hadrons. 

For the parton to lose a finite fraction of its forward 
energy, multiple gluon emissions are required. Each such 
emission at a given opacity is assumed to be independent 
and a probabilistic scheme is set up, wherein, the jet loses 
an energy fraction AE in n tries with a Poisson distribu- 
tion 22,47], 



-<N B ) 



11 i=l 



(10) 



where, (N g ) is the mean number of gluons radiated per 
coherent interaction set. Summing over n gives the proba- 
bility P(e) for an incident jet to lose a momentum fraction 
e due to its passage through the medium. 

This probability distribution is then used to model a 
medium modified fragmentation function, by shifting the 



energy fraction available to produce a hadron as well as 
accounting for the phase space available after energy loss 
(The fragmentation function used is a vacuum fragmenta- 
tion function) . The medium modified fragmentation func- 
tion is thus defined as 22,(57], 



f i D 

D(z,Q 2 )= / deP(e) — 
Jo 



1 - e 



(11) 



The above, modified fragmentation function is then used 
in a factorized formalism as in Eq. (fT]) to calculate the 
final hadronic spectrum. Additional details related to the 
computation of the energy loss probability distribution are 
given in Ref. [34] . 



4.3 Finite temperature field theory scheme 

In subsection 13. 3( the computation of the rates of parton 
splitting and merging in a thermalized deconfined medium 
were calculated. Based on these rates, a Fokker-Plank 
equation was motivated which computed the change in 
the distribution of hard partons with time spent prop- 
agating through a medium. If the initial distribution is 
taken from the cross-section to produce a hard parton as 
in Eq. (fT]) and the time spent in the medium is estimated 
based on the production point and the direction of prop- 
agation, this equation will yield the distribution of hard 
partons as they exit the deconfined medium. In the AMY 
scheme, these partons are no longer expected to interact 
with the hadronic plasma and thus do not lose energy in 
the hadronic phase. 

The final hadron spectrum at high px is obtained by 
the fragmentation of jets in the vacuum after their passing 
through the medium. In this approach, one calculates the 
medium modified fragmentation function by convoluting 
the vacuum fragmentation functions with the hard par- 
ton distributions, at exit, to produce the final hadronic 
spectrum 28 , 



Dj(z, r±, i 



dp r ^D^(z')P(p f \ Pj ,r ± ,4>). (12) 



In the equation above, the sum over j' is the sum over 
all parton species. The two momentum fractions are z = 
Ph/Pj and z' — Ph/Pj', where pj and pji are the momenta 
of the hard partons immediately after the hard scatter- 
ing and prior to exit from the medium and ph is the final 
hadron momentum. The quantity P(j>y \pj, r±, ip) repre- 
sents the solution to Eq. ([5]), which is the probability of 
obtaining a given parton j' with momentum py when the 
initial condition is a parton j with momentum pj. The 
above integral depends implicitly on the path taken by 
the parton and the medium profile along that path, which 
in turn depends on the location of the origin r± of the jet 
and its propagation angle <f> with respect to the reaction 
plane. Therefore, one must convolve the above expression 
over all transverse positions r± and directions </>. Details 
of this procedure are presented in the companion paper of 
Ref. [31. 
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5 Discussions and conclusions 



6 Acknowledgments 



In these proceedings, the different underlying theoreti- 
cal mechanisms used in some of the prevalent jet energy 
loss calculations have been outlined. Specific attention was 
paid to how each jet resolves the medium and on the prop- 
erty of the medium which controls the modification of the 
hard jet. In all cases, the jet modification formalisms may 
be reduced to a form which depends on only one parame- 
ter: this is the transport coefficient q defined as the trans- 
verse momentum gained by a hard parton per unit length 
traversed in a dense medium. While in the ASW formal- 
ism, q is the sole tunable parameter, in the HT formalism, 
q depends on the gluon field strength correlation [Eq. |4j 
and may be calculated from a knowledge of the tempera- 
ture and the coupling constant a s in the AMY formalism. 

In all cases, the jet is assumed to fragment outside 
the medium. As a result, all formalisms use a medium 
modified fragmentation function, which uses a vacuum 
fragmentation function as input. While in the ASW and 
the HT formalisms, the modification is computed in both 
deconfined and confined phases, due to the assumptions 
made in the AMY formalism, the modification in this for- 
malism occurs only in the deconfined phase. The modifi- 
cation in the confined phase is assumed to be small and 
ignored. While both the HT and the ASW formalisms in- 
clude contributions from interference with vacuum radia- 
tion, these are ignored in the AMY scheme. The AMY ap- 
proach however includes contributions from thermal feed- 
back which has so far not been straightforwardly included 
in the HT and ASW formalisms. The consistent setup of 
the AMY formalism also allows for the most natural exten- 
sion to include elastic energy loss [41] . While in the strict 
interpretation of heavy scattering centers in the ASW (and 
GLV) formalism, elastic energy loss is identically zero, the 
inclusion of elastic loss requires additional assumptions 
about the medium in the HT approach. Extensions to the 
GLV formalism to include mobile scattering centers, and 
thus, both include elastic energy loss [33] and modify the 
formulation of radiative energy loss 48J are currently un- 
derway. Similar extensions in the HT approach are also be- 
ing carried out |49j . However, given the incomplete setup 
in different formalisms, elastic energy loss was not dis- 
cussed in these proceedings, nor will be included in the 
realistic comparisons presented in Ref. [34] , 

While the description of the different formalisms in 
these proceedings have not included the effect of a dy- 
namical medium, realistic calculations of jet modification 
in heavy-ion collisions do include such effects. The modifi- 
cation depends on the path traversed by a given jet. This 
in turn depends on the origin of the jet and the direction of 
travel in the medium. The details related to this problem, 
as it applies to the different formalisms will be presented 
in the companion paper [34) . As a result, comparisons to 
experimental data will also be carried out in this reference 
as well. 
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